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Abstract. We determine the fundamental stellar parameters of altogether 20 double-mode RR Lyrae (RRd) stars from the 
Galactic Field and the Large Magellanic Cloud, following our approach employed on the field RRd star BS Comae 1 1]. The 
stars were selected to cover wide ranges of periods and period ratios, implying diverse stellar parameters. From the possible 
observed quantities we use only the periods and determine stellar parameter combinations that satisfy both current helium- 
burning evolutionary models and grids of linear non-adiabatic purely radiative pulsational models. Thus, the periods of an 
object determine a sequence of solutions for its mass, luminosity, effective temperature and metallicity, parametrized by the 
time elapsed from the zero age horizontal branch. The derived sets of solutions yield various important theoretical relations 
between the physical parameters of the stars which are, in the case of some parameter combinations, nearly independent of 
the age. We get very tight simple linear relations between logPo and log/?, logp, logg and the Wesenheit index W{B — V). 
This latter period-luminosity-color relation is in good agreement with the one derived on empirical basis [2j and calibrated by 
Baade-Wesselink results |3]. 
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INTRODUCTION 

Double-mode RR Lyrae (RRd) stars pulsate simultane- 
ously and stably in the first two radial modes with a pe- 
riod ratio close to 0.75. Their two precisely observable 
periods make them a particularly important subtype of 
RR Lyrae stars by supplying an additional constraint on 
their structure. Eventually, the double-mode nature of- 
fers the possibility of determining the fundamental pa- 
rameters of these stars from radial pulsation models. The 
solid theory of linear radial pulsation has long been used 
for this purpose by fitting models simultaneously to the 
periods of both pulsation modes. By measuring or giving 
an estimate on the chemical composition and effective 
temperature of the star we can derive its mass and lumi- 
nosity. The consistency of these values can be confirmed 
by independently obtained distances [see, e.g., 4]. 

Nonlinear hydrodynamical modeling of double-mode 
pulsation has also became possible in the past decade Q 
1^. Sustained double-mode pulsation can be reproduced 
in a consistent range of physical parameters, a star-by- 
star modeling, however, is still not available. 

For a fully consistent RRd model it is important to 
take horizontal branch (HB) stellar evolution also into 
account. In a novel study of this endeavor, to involve evo- 
lutionary constraints in explaining the observed period 
distribution of RRd stars of the Large Magellanic Cloud 
(LMC) by means of physical parameter ranges, |7] calcu- 
lated pulsation model sequences along theoretical tracks 
of HB evolution to compare the theoretical period dis- 
tribution with the observed one. In the present study we 



follow a different approach by taking the observed peri- 
ods as an input to our calculations and search for a set 
of physical parameters that satisfy the combined param- 
eter space of pulsational and evolutionary models. Thus 
obtaining strongly constrained values for the fundamen- 
tal parameters of a representative sample of RRd stars, 
we can probe the full parameter ranges they populate, as 
well as use them to deduce theoretical relations among 
their physical parameters. 



METHOD AND DATA 

Our goal is to determine the fundamental parameters of 
an RRd star solely from its periods by restricting the pa- 
rameter space of pulsation models by the constraints of 
evolutionary models. The basic assumption behind this 
approach is that there is a good level of consistency be- 
tween the two kinds of models. We obtain the fundamen- 
tal physical parameters of a star by the same procedure as 
described in 1 1], and hereby we only briefly summarize 
the essentials of the method. 

On the one hand we calculate a grid of linear pulsa- 
tion models fitting the observed pair of periods. These 
models give stellar mass and luminosity for specified ef- 
fective temperature and metallicity values. On the other 
hand, isochrones of Helium-burning models also deter- 
mine mass and luminosity as a function of effective tem- 
perature, metal content, and age. To obtain the set of pa- 
rameters that satisfy both models we have to calculate the 
intersection of these two restricted parameter-subspaces 
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FIGURE 1. An example for the combined pulsational and 
evolutionary solution sequences calculated for an RRd star in 
the LMC. Age dependence of the various physical parameters 
are shown up to log/ = 7.86 (~ 72.5Myr). 
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FIGURE 2. Petersen diagram showing RRd stars in the 
LMC, taken from the MACHO catalog I?, dots], together with 
all known Galactic Field RRd stars taken from the literature 
[3,111111,111113111113111 crosses]. Stars analyzed is this 
study are denoted by circles. 



(formally vector-vector functions that we sample on a 
grid). By fixing the periods to the observed ones, evo- 
lutionary models have one variable in surplus, the age, 
ie. the time elapsed from the zero age horizontal branch 
(ZAHB), which will be a free parameter of the solution. 

We used the canonical HE evolution models of ||8t] 
with solar-scaled heavy element mixture. These mod- 
els have a fine mass and time spacing, suitable for the 
detailed investigation of the instability strip. We ob- 
tained isochrones by interpolating the tracks to a uniform 
mass and time base. We generated linear, non-adiabatic 
(LNA), fully radiative pulsation models. The pulsation 
model grid was calculated for the same chemical compo- 
sitions as those of the evolutionary tracks, and with a fine 
temperature step of 50 K. The deep and dense sampling 
of the stellar envelopes resulted in an accuracy of ~ 10^^ 
in the period fit. See JH] for further details on the model 
construction. 

As a result of the method, for a pair of observed peri- 
ods we obtain a 1-D solution sequence giving the effec- 
tive temperature, metallicity, stellar mass, and luminosity 
as a function of the time parameter Thus, each point of a 
sequence determine combinations of physical parameters 
corresponding to different stages of evolutionary tracks 
of different masses and metallicities. We also calculate 
synthetic colors from the parameter sequences using the 
stellar atmosphere models of |9]. Fig.[T]gives an example 
solution for the various physical parameters. We note that 
in the absence of any additional reliable observational 
constraint, the models themselves alone do not pose any 
limit on the time parameter. However, we can still put 
rough boundaries on the solution sequences by keep- 



ing them within reasonable temperature ranges or omit- 
ting more advanced evolutionary stages that are more 
unlikely to be observed due to fast parameter changes. 
(Also, solutions that would require extremely low metal- 
licities outside the evolutionary model grid can techni- 
cally confine the possible ages.) Thus, in most cases we 
stop scanning in age between logf = 7.85 and 7.9. 

Figure |2] shows the Petersen diagram of all known 
Galactic Field RRd stars and those found in the LMC 
by the MACHO project. The ranges of the Pi/Pq vs. 
Po distribution populated by these stars encompass the 
vast majority of all known double-mode RR Lyrae stars 
in any stellar system. We selected a sub-sample of 20 
stars (3 from the Field: AQ Leo, BS Com, and CU Com, 
and 17 more from the LMC) to survey full ranges of 
physical parameters, as implied by the periods and period 
ratios. Sequences for fundamental physical parameters 
were obtained for each of the stars in the way discussed 
above. We present the analysis of the derived sets of 
solutions in the following section. 

RESULTS 

Based on the Petersen diagram, the physical parameters 
obtained for the 20 objects investigated in this work are 
representative for all known RRd stars (and their lumi- 
nosities, presumably, for a dominant fraction of single- 
mode RR Lyrae stars). Therefore, we use them to mea- 
sure the fundamental parameter ranges populated by 
these objects and probe important theoretical relations 
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FIGURE 3. The loci of the RRd stars analyzed in this paper 
on the HRD. Solution sequences with free time parameter are 
shown as thick lines. Values corresponding to zero ages are 
plotted by circles. Helium-burning evolutionary model tracks 
of 1 8] are shown by thin gray lines for [Fe/H] values of —1.3, 
— 1.5, and —1.8 (redward parts of tracks for masses lower than 
0.64, 0.66, and 0.70, respectively, are not plotted for more 
perspicuity). 



defined by the solutions of our method. 

Figure |3] shows the distribution of the solutions on 
the Herzschprung-Russell diagram (HRD) together with 
various evolutionary tracks for reference. The expected 
general trend, that stars with lower metallicity are cooler 
and brighter, is clearly exhibited. The gradual topological 
change of the sequences originates from the smoothly 
changing sizes and shapes of the blueward loops in the 
constituting tracks, mainly as a function of metallicity. 
As it is also shown in Fig.lT] the temperature dependence 
of a solution increases considerably for higher values of 
age (i.e., after reaching the high temperature extrema - 
the 'blue noses' - of the constituting tracks). Therefore 
we would need further specific information on each star 
(i.e., on their spectral energy distributions) to put more 
significant constraints on their temperatures. However, 
according to the combined models, RRd stars populate 
an overall temperature range of at least 300 K, as it is 
clear from Fig.|3] 

Figure |4] shows the solutions on the B — V vs. V — I 
two-color diagram. Stars are concentrated into a very 
narrow overall range of color index pairs. Here we note 
that by a very precisely measured pair of color indices 
one would be able to further constrain the time parame- 
ter, thus reducing the ranges of uncertainty in all other 
parameters, as it was successfully done in the case of 
BSCom [see 1]. 

In Fig. |5] we examine the luminosity and mass as 
functions of the metallicity. Values corresponding to the 
ZAHB states show a tight and, apparently, somewhat 
non-linear correlation with [Fe/H] in both cases. Allow- 
ing later stages of HE evolution complicates the topol- 
ogy, but the models set minimal luminosity and maxi- 
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FIGURE 4. Two-color representation of the RRd solution 
sequences. Notation is the same as in Fig. [3] 



mal mass values at any specified metal content. Accord- 
ingly, for specified mass the models give a minimal lu- 
minosity, and, vice-verse, for a fixed luminosity, they put 
a constraint on the maximum value of the mass. Since 
we selected our stars to cover nearly the full RRd period 
ranges of the LMC (and by this selection, that of almost 
all known RRd stars as well), the solutions inherently 
yield the full theoretical [Fe/H], mass, and My ranges 
spanned by the constituting objects. Although the indi- 
vidual sequences suffer from an uncertainty of roughly 
0.1 — 0.3dex in metallicity (depending on the allowed 
range of the time parameter), the upper limit for the over- 
all [Fe/H] range turns out to be ^ 0.9 dex with extrema 
of approximately —1.3 and —2.2. We note that, as im- 
plied by the Petersen diagram (ie., the short period end is 
more populated in the case of the LMC), the majority of 
the stars are in the upper part of this metallicity range. 
By direct spectrophotometric measurements of 6 RRd 
stars in the LMC, [18] found a wider metallicity range 
of (-2.28,-1.09). Besides the individual errors of the 
few measurements and the sensitivity of the AS method 
to different calibrations, this moderate difference may re- 
sult from the possible non-solar heavy element mixture 
of these stars (we recall that throughout this paper solar- 
scaled mixtures of heavy elements are used). 

The stellar masses corresponding to our metallicity 
domain span in a large O.IVM© range from 0.64 M© 
to 0.81 Mq. As for absolute magnitudes, the luminosity 
interval permitted by the solutions on (see Fig. |5]l can 
be translated into an My range of (0.46,0.70) which, as 
we will see later, is in a good agreement with Baade- 
Wesselink results. 

We also computed \ogg, log 7?, and log p values for all 
solution sequences. These parameters are highly degen- 
erate, i.e., very insensitive to the time parameter. They 
also show very tight, linear correlations with the funda- 
mental period. Figure |6] shows these simple theoretical 
relations provided by our results. It is remarkable how 
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FIGURE 5. Luminosity-metallicity and mass-metallicity plots of the RRd stars studied in tliis paper. For notation see Fig. [3] 



tiny individual parameter ranges are allowed by the so- 
lutions. Thus, for the sake of simplicity, linear formulae 
given on the same figure were obtained by fitting to the 
zero age solutions only. 

One of the tightest empirical relations for RR Lyrae 
stars is the period-luminosity-color (PLC) relation, or in 
other words the period - reddening-free Wesenheit-index 
relation. We can directly compare it with the one derived 
from the synthetic My and color index values of RRd 
stars calculated above. Figure |7] shows the theoretical 
PLC relation for B and V colors stretched by our RRd 
model results. In this representation of the solutions we 
notice a vast insensitivity again to evolution effects. For 
comparison we show the empirical relation derived from 
a large number of RRab stars from Galactic globular 
clusters (GCs) by |2], together with the 3 <J range of the 
calibration sample. The zero-point of this relation was set 
by the Baade-Wesselink (BW) results for Galactic Field 
RRab stars (also shown in Fig. 13 as given by |3]. The 
zero points of the two relations show a fine agreement. 
Although the RRd sample clearly defines a significantly 
steeper slope, it does not exceed the 3 (7 ranges of the GC 
sample. We note that the situation is very similar in the 
case of the V and / colors. 

The reason for the slope difference is not yet clear 
and may be caused by a combination of different effects. 
Since the PLC relation is governed by two parameters 
(i.e., metallicity and effective temperature) out of four of 
the pulsation equation, it has an intrinsic scatter within 
which the subset of RRd stars could follow their own, 
tighter relation, determined by the specific stellar param- 
eter subspace they populate. 

We can also speculate on the following. In the deriva- 
tion of the empirical PLC relations by |01, no direct 
calibration (i.e., using independently measured absolute 
magnitudes) was made. Instead, the relative dereddened 
distance moduli of the constituting globular clusters were 



obtained as output parameters of the least-squares mini- 
mization procedure, together with the regression coeffi- 
cients. Thus, a uniform slope was obtained for all stars, 
nicely fitting the subsamples of the individual clusters, 
and the zero-point differences were interpreted as the re- 
sult of the relative distances. Later [19] found a some- 
what steeper Wiyj) - logPo relation for LMC RRab 
stars than that of but the scatter was much higher. 
The HE stars of a GC have generally little if any metal- 
licity dispersion. Therefore, if an additional slight, intrin- 
sic metallicity dependence of the Wesenheit index does 
rather manifest in zero point than in slope differences 
between subsamples of stars with constant metallicity, 
its effect could merge with that of the relative distances 
and remained hidden in the above case. Furthermore, if 
this is so, then our RRd sample will yield a steeper rela- 
tion because it scans a considerable, nearly 1 dex range in 
[Fe/H], which, in this case, is more correlated with logPo- 
However, we would need very precise independent rela- 
tive distance determinations for the calibrating clusters 
to give a more secure assessment of this problem. 

CONCLUSIONS 

We presented a method and its applications for stellar pa- 
rameter determination of double-mode RR Lyrae stars. 
The method combines theoretical model results of lin- 
ear stellar pulsation and canonical HB stellar evolution. 
Based on the periods of an object, the joint constraints of 
the models allow sequences of fundamental stellar pa- 
rameters with age as a free parameter. Additional ob- 
servables like color index or metallicity can be used to 
put further constraints on the solutions. By using a repre- 
sentative sample of objects to probe the observed period 
ranges of the known RRd's with our method, we obtained 
the theoretical ranges of physical parameters populated 
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FIGURE 6. Surface gravitational acceleration, stellar radius, and stellar density as a function of the period of the fundamental 
mode. Notation is the same as in Fig.jS] Lines show linear fits to the zero age solutions, yielding the formulae given in each panel. 
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FIGURE 7. PLC relation for the RRd sample of this study 
are shown with the same notation as in Fig. [3] Crosses denote 
BW results of |3]. Solid line shows the linear fit to the RRd 
data. PLC relation of |2] with zero-point fixed by the BW 
sample is delineated by a dashed line, with a shaded area 
illustrating the 3 o ranges of its calibrating sample. 



by these objects. We also extracted important theoreti- 
cal relations from the derived sets of stellar parameters 
among which we found nearly age-independent ones. We 
derived a very tight linear PLC relation with a magnitude 
zero-point being in a fine agreement with that of Baade- 
Wesselink results. We found a slope difference, which 
might be explained by a hitherto hidden metallicity de- 
pendence of the relation. 
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